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A NEW RESULT ON THE ORIGIN OF THE 
EXTRAGALACTIC GAMMA-RAY BACKGROUND 



ABSTRACT 

In the paper, we continually use the method of image stacking to study the 
origin of the extragalactic gamma-ray background (EGB) at GeV bands, and 
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find that the Faint Images of the Radio Sky at Twenty centimeters (FIRST) 
sources undetected by the Large Area Telescope on the Fermi Gamma-ray Space 
£3 ■ Telescope can contribute about (56±6) % of the EGB. Because the FIRST is a 

flux limited sample of radio sources with incompleteness at the faint limit, we 
consider that the point-sources, including blazars, non-blazar AGNs, starburst 
galaxies, could produce a much larger fraction of the EGB. 

in 

iy-j ■ Subject headings: gamma rays: diffuse background — methods: statistical — quasars: 

CO ■ general — galaxies: starburst 

O 
CO 

1 

1. Introduction 

The extragalactic gamma-ray background ( EGB) at GeV band s was first detected by 



the satellite Small Astronomy Satellite- 2 (SAS-2; iFichtel et al.lll975l ). and its s pectrum was 



measu red with good accuracy by Fermi (also called isotropic diffuse background; lAbdo et al. 



2010|). Based on the first-year Fermi data, the EGB has been found to be consistent with a 



featureless power law with a photon index of ~2.4 in the 0.2-100 GeV energy range and an 
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integrated flux (E>100MeV) of 1.03xl0~ 5 ph oton cm 
E>lGeV is 4xl0~ 7 photon cm -2 s _1 sr _1 (see 



-2 — 1 

s sr 



Abdo et al.ll2010d ). 



The integrated flux with 



The origi n of the EGB is one of the fundamental unsolved problems in astrophysics (see 
Kneiskell2008l . for a review). The EGB could originate from either truly diffuse processes or 
from unresolved point sources. Tru ly diffuse emission can arise from s o me processes such 
as th e annihilation of dark matter (jAhn et al.l 120071 : ICuoco et al.l l2010l ; iBelikov fc Hooper 
20101 ). the emission of high energy particles acc elerated by intergalact ic shocks which are 
produced during large scale structure formation (IGabici fc Blasii .1120031 ) etc. 



Blazars, including BL Lac objects, flat spectrum radio quasars, or unidentified flat 
spectrum radio sources, represent the most numerous population detected by the Energetic 
Gamma Ray Experim ent Telescope (EGRET) on the Comyton Gamma Ray Observator y 
(IHartman et al.lll999h and by the Fermi-L AT (lAbdo et alihoiObl : lAckermann et alihoill ). 
Therefore, blazars undetected by the EGRET or Fermi-LAT are the most likely candidates 
for the emission of the EGB. Many authors have studied the luminosity function of blazars 
and showed that the contribution of undetected blazars to the E GB could be in the range from 



20 % to 100 % flStecker & Sa 



2008 



2011 



Kneiske fc Mannheim 



Stecker fc Ventersll201ll ). However 



anion 



2008; 



1996; lNarumoto fe Totani 20061; Dermer 



Inoue fc Totani 



Abdo et al. 



2007 



Cao fc Bai 



2 0o3lLi&; Caol 120111 

(j2010d ) built a source count distribu- 



Malyshev h Hogg 



tion at GeV bands and yielded that blazars undetected by the LAT can contribute about 
23 % of the EGB. They ruled out blazars producing a large fraction of the EGB. 

Non- blazar radio loud active galactic nucle i (AGNs) can also contribute a fraction to 
the EGB ( Bhattacharva et al. 12009c llnoud 120111 ). since some of them have been found to be 



7-ray sources ( lAbdo et al. 



2010dl ) and they outnumber blazars. 



Starburst galaxies, which have intense star-formation, are expected to have high super- 
nova (SN) rates. SN remnants are believed to accelerate primary cosmic rays (CRs) protons 
and electrons. The high SN rates in starburst galaxies imply high CR emissivity. When 
high energy CR protons collide with interstellar medium (ISM) nucleons, they create pions 



decayin g into secondary electrons and positrons, 7-rays, and neutrinos. lAckermann et al. 
()2012aJ) have found quasi-linear scaling relations between 7-ray luminosity and both radio 
continuum luminosity and total infr ared luminosity. To data, several starburst galaxies have 



been d etected by the Fermi-LAT (lAbdo et al.l l2010at iNolan et al.l 120121 ; lAckermann et al. 



2007 



2012a ). Therefore, starbursts can also contribute a fraction of the EGB (|Thompson et al. 



Bhattacharva fc Sreekumarl 120091 ; lLacki et al.l 1201 ll ; iMakiya et al.l l201ll ) . 



Nevertheless, normal galaxies (Bhatt acharva fc Sreekumarl 120091 ; Fields et al.ll2010l ) or 
radio-quiet AGNs ( llnoue et al.ll2008l ; llnoue fc Totanill2009l ) can also contribute a fraction of 
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the EGB. 



We introduced a new method of image s tacking to directly study the undetected but 

20111). A pplying the method to the Australia Tele- 



possible 7-ray point sources (IZhou et al. 



scope 20 GHz Survey fAT20G: lMurphv et al. 



20101 ) sources undetected by the LAT, we found 



that these sources contribute about 10.5 % of the EGB in the 1-3 GeV energy ranges. 

In this paper, we continually use the method to study the origin of the EGB, and find 
that the FIRST sources undetected by the LAT can contribute about (56±6) % of the EGB. 
We also discuss the implications. 



2. FIRST 



The Faint Images of the Radio Sky at Twenty centimeters (FIRST; iBecker et al.lll995l ) 
program is using the Very Large Array (VLA) to produce a map of the 20 cm (1.4 GHz) sky 
with a beam size of 5. "4 and an rms sensitivity of about 0.15 mJy beam -1 . The accuracy of 
source position is better than 1". The survey has covered an area of about 9,055 deg 2 in the 
north Galactic cap and a smaller area along the celestial equator, corresponding to the sky 
regions observed by the Sloan Digital Sky Survey (SDSS). At the 1 mJy source detection 
threshold, the source surface density is ~ 90 deg -2 , and the catalo g includes ~ 8 x 10 5 
sourced_|, in which about 30% of them have counterparts in the SDSS (llvezic et al.l 120021 ). 



The FIRST catalog lists two types of 20 cm continuum flux density, e.g. the peak 
value, -F peak , and the integrated flux density, F int . These measurements are derived from 
two-dimensional Gaussian fitting for each source, where the source maps are generated from 
the co-added images from 12 pointings. 



Ivezic et al.l (120021 ) found that about 30 % of FIRST sources was positional association 



within l."5 with an SDSS source in 1230 deg 2 of sky. The majority (83%) of the FIRST 
sources identified with an SDSS source brighter than r* = 21 are optically resolved. About 
30 % of them are radio galaxies in which emission-line ratios indicate AGNs; the others 
are starburst galaxies. Nearly all optically unresolved radio sources have nonstellar colors, 
indicating quasars. Because there is no significant difference in the radio properties between 
FIRST sources with and without optical identifications, the majority of unmatched FIRST 
sources could be too optically faint to be detected in SDSS images, and the fractions of 
quasars and galaxies are roughly the same for the two subsamples. 



1 It is available online http://sundog.stsci.edu/nrst/catalogs/readme_08jull6.html 
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Because two subsamples of the FIRST catalog are 7-ray emitters, this catalog maybe a 
5000! tracer for the undetected 7-ray point sources. 



Method 



For a sample of poss ible 7-ray point so urces undetected by the Fermi due to their 
faint fluxes, soft spectra (jAbdo et al.l |2010eJ) or the source confusion (IStecker h Venters 
2011). we can stack a large number of them to improve the statistics ( IWhite et al.l 12007 ; 
Ando fc Kusenkd 120101 ). We have use d the Maximum Likelihood (ML) method to derive 
the fluxes of the stacked point sources (IZhou et al.l 1201 ll ). But for a very large sample, this 
method is very time consuming. Therefore, in this paper, we use a rough but simple ML 
method to derive the fluxes. 

The photongj we used in our analysis are taken during the period of 2008 August 4 
(15:43 UTC) - 2011 October 20 (23:33 UTC), about three years. During most of this time, 
Fermi was operated in sky-scanning survey mode (viewing direction rocking north and south 
of the zenith on alternate orbits). Time intervals flagged as 'bad' (a very small fraction) and 
the period of the rocking angle of the LAT greater than 52° was excluded. Only the photons 
in the 1-100 GeV e nergy range with small 68 % contain ment radius (better than 1°) and 



little confusion (see lAtwood et al.ll2009l : lAbdo et al.ll2009l ) are used. These photons are also 



needed to satisfy the standard low-background event selection (termed "Source" class events) 
corresponding to the P7V6 instrument response functions in the present analysis. The effect 
of Earth albedo backgrounds was greatly reduced by removing photons coming from zenith 
angles < 100°. In this procedure, the tools of gtselect and gtmktimeU are used. 

Stacking the images of the sources, we collect all photons that are at most 1° away from 
any source of our sample and then record their angular distance (6i, in units of deg) between 
the photon and the source. The 1° is enough despite some photons attributed to the central 
stacked point sources are not in the stacked image. Due to the faintness of the stacked point 
sources, in the outer of the stacked image, the signal of the stacked point sources is very 
weak. 

In order to minimize the influ ence of strong poin t sources, the sources away from any 
second Fermi-LAT catalog (2FGL; iNolan et al.ll2012l ) sources less 2° (3° for 2GFL sources 



2 http://fermi. gsfc.nasa.gov/cgi-bin/ssc/LAT/LATDataQuery.cgi 

3 Thcse and other tools we used in next are parts of most recent Fermi-LAT Science Tools, version v9r23pl, 
which are accessible at |http://fermi.gsfc.nasa.gov/ssc/data/analysis/scitools/overview.html| 
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with fluxes larger than 10~ 9 photon cm -2 s _1 in the 1-100 GeV energy range) are not used. 
In order to minimize the influence of strong emission from the Galactic plane, only the 
sources locating at high Galactic latitudes, e.g. \b\ > 15°, are used. 

There are a s mall number of spur ious sources that are sidelobes of nearby bright sources. 



In recent catalog. iBecker et al.l (120031 ) uesd P(S) to indicate the probability of that the source 
is spurious (most commonly because it is a sidelobe of a nearby bright source). The mean 
value of P(S) is 0.088, which indicates th at about 8.8 % of the sour ces are sidelobes. They are 



concentrated around bright sources (see lWhite et al.l 119971 120131 for detail), and should be 



removed first. In Figure 1 we show the number density of all FIRST sources and the sources 
having P(S) < 0.1 around the FIRST sources with F pea k >100 mJy, respectively. A strict 
P(S) cutting, such as P(S)<0.1, is not appropriate. We randomly remove some sources with 
the probability of their P(S) except the sources with P(S)=0.014 (e.g., the minimum value 
of P(S), for about 72% of all sources). This procedure does not try to remove sidelobes as 
much as possible, but to make the distribution of the rest sources around the bright sources 
as uniform as possible. The number density of the rest sources around the FIRST sources 
with -Fp ea k >100 mJy is also shown in Figure 1. The distribution is still not uniform enough, 
we will investigate the effect of uniform distribution in the following section. 

Many sources in the FIRST catalog are not independent o bjects but are components 



of a single object with complex morphology (IWhite et al.lll997l ). In Figure 2, we show the 



mean source number distribution of other FIRST sources around one FIRST source. It is 
found that the possibility of finding another FIRST source near a FIRST source, especially 
closer than ~ 50", is larger than one expected from a random source distribution. Therefore, 
when a group sources fall within 50" of their nearest neighbors, we treat them as a single 
object. Nearly 30 % of the FIRST sources belongs to such sources. The cutoff of 50" is a 
trade-off between the completeness and contamination of the sample. For a larger cutoff, 
more complex sources are merged, but more independent sources are also merged. The use 
of 50" cutoff will produce small uncertainty introduced by the misclassified sources, and will 
be discussed in the following section. 

The photon number density profile of the stacked image of FIRST sources with -F pca k > 1 m Jy 
is shown in Figure 3. It is shown that the center of the stacked image has higher photon 
number density and the stacked point source has nonzero flux. 

We use a e (6) to further prove the presence of the point source, which is the confidence 
level of stacked image representing photons within the radius of 9 more than the expectation 
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in the uniform density. It is defined as 



n(9) - §e 2 



N_m 

e 2 ° 



where N(9) is the number of photons with 9i < 9, N is the number of photons with 9i < G, 
here = 1°. The a e (9) is presented in Figure 4. The excess of photons is very obvious 
within 9 ~ 0.3°. 

we apply the maximum likelihood (ML) method to derive the flux of the stacked point 
source. The likelihood is the probability of the observed data for a specific model. For our 
case, it is defined as 

33 

L = HlH, (2) 



i=\ 



where 



1 



Pi 



(Ni - M % ) 



V^N^i ^ J (3) 

is the normal probability of observing Mi counts in i-th bin when the number of counts 
predicted by model is JVj. The logarithm of the likelihood is conveniently calculated: 



33 



i„L = yj 



i=l 



-0.51n(27r) -0.5 In N t 



(Ni - M t ] 
IN 



(4) 



For simplicity, in our model there are only two components, e.g., the central stacked point 
source and diffuse background source, and 



Ni = N s n si + N b n bi , 
N = N S + N b , 



(5) 
(6) 



where N s and N b are the photon numbers attributed to the central stacked point source and 
the diffuse background source, respectively. n S i and nu are the normalized photon number 
distributions of two class sources. 

In order to determine nu, we stacked 10 8 imaginary sources which is isotropically dis- 
tributed on the sky. In fact, only ~ 4 x 10 7 sources far away strong sources are ultimately 
stacked. Because the sources we stacked here are not true 7-ray point sources, the stacked 
image only contains diffuse backgrounds. The photon number density profile of the stacked 
image is shown in Figure 5. It is found that the diffuse background is not uniform and has 
higher density in the outer of the stacked image, implying that the profile is affected by 
the 2FGL sources near the stacked image. In Figure 6 we show the photon number density 
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profile of the stacked image which is extended to 2°, indicating that the effect of the nearby 
point sources is much obvious. 



s 1 in 



In order to determine n si , we simulate a source with a flux of 5 x 10 6 photon cm 2 
the 1-100 GeV energy range and the coordinate of (RA, DEC) = (190°, 30°fl It has 
a power law spectrum and the photon index is 2.4, which is a typical value of the detected 
sources. The tool gtobssim has been used in this procedure. The photon number density 
profile of this source is shown in Figure 8. 

We use the likelihood ratio to test the hypothesis. The point-source "test statistic" (TS) 
is defined as 

TS = -2(lnL -hiL 1 ) (7) 

where L\ and L are the likelihood with and without point source. The TS of each source 
is related to the probability of that such an excess is obtained from background fluc- 
tuations alone. The probability distribution in such a situation is not known precisely 



(jProtassov et al.ll2002l ). However, we only consider positive fluctuations, i n which each fit 



ting involves one degree of freedom, and according to Wilks's theorem (jWilksl Il938l ) the 
probability with at least TS is close to that of the x 2 /2 distribution with one degree of 



freedom. Therefore the detected significance of a point source is approximately yTSa (see 



Mattox et alJll99q ). 



4. Result and Discussion 

For the stacked point source of the FIRST sources with F pea k >1 mJy, the estimated N s 
is 2.51 x 10 4 , in which its TS is 43.0, corresponding to a significance of ~ 6.7cr. The x 2 over 
degree of freedom (dof) is x 2 /dof = 30.7/31, indicating that the model gives a satisfactory 
fit to the data. 

N s is proportional to the flux and exposure of the stacked point source. The exposure 
is the mean of the exposures of all stacked sources. But the calculation of the exposures^ 
for all stacked sources is very time consuming, and the exposure of all sky is relatively 
uniform, owing to the large field-of-view and the rocking-scanning pattern of the sky survey 



(INolan et al.ll2012f ). we use the mean exposure of about 1 % randomly selected FIRST sources 
to represent them. The exposures of the stacked point source and the simulated source 
are shown in Figure 7, which are roughly the same. The simulated source has a flux of 



This point has similar exposure to the stacked point sources, see Figure 7. 
5 Calculated using the tool gtpsf in Fermi-L AT Science Tools. 



5 x 10~ 6 photon cm" 2 s _1 and ~ 4.1 x 10 5 photons within 1°. The stacked point source has 
~ 2.51 x 10 4 photons within 1°, corresponding to a flux of ~ 3.09 x 10~ 7 photon cm -2 s _1 . 
The FIRST covers an area of about 9055 deg 2 , in which 321556 out of 640919 sources are 
stacked, indicating that the stacked point sources cover an area ( i.e. the area covered by 
the FIRST sample excluding all holes produced in the Fermi map by excluding photons in 
the vicinity of 2FGL sources.) of about 4543 deg 2 . In this sky region the total flux of the 
EGB is ~ 5.7 x 10~ 7 photon cm -2 s _1 , implying that these sources can contribute ~ 56 % 
of the EGB. 

A decrease of 0.5 from i ts maximum value i n InL corresponds to 68% confidence (1 a) 



region for the parameter (see iMattox et al.lll996l ). We use this variance to estimate the error 



of the flux and find that 1 a error is 0.33 x 10 7 photon cm 2 s 1 , implying that the FIRST 
sources contribute about (56± 6) % of the EGB. 

The effective area is an important factor of uncertainties. The current estimate of the 
remaining systematic uncertainty is 10 % a t 100 MeV, decreasing to 5 % at 560 MeV and 
increasing to 10% at 10 GeV a nd above (lAckermann et al.l l2012bl ). This uncertainty is 



uniformly applied to all sources ( Nolan et al.l 120121 ) . The error of the fraction of FIRST 



sources contributed to the EGB is much smaller. 

In order to verify the effectiveness and accuracy of our method, we do the Monte Carlo 
simulation using the tool gtobssim. The simulating time is ~ 10 8 s, equaling the time of real 
data we used. We simulate 7.8 x 10 5 point sources with fluxes of 10~ 12 photon cm -2 s _1 in 
the 1-100 GeV energy range. The flux distribution does not affect our results, which only 
depend on the total flux. They have power law spectra, in which the photon index of each 
source is drawn from a Gaussian distribution with a mean of 2.40 and 1 a width of 0.24, 



which is the same with the intrinsic distribution of photon indices (jAbdo et al.ll2010d ). The 
coordinates of each source are randomly drawn to produce an isotropic distribution on the 
sky between right ascensions of 95° and 270°, declinations of -10° and +70°, which covers the 
main area observed by the FIRST survey. The simulated sources have similar source surface 
density with the FIRST sources, in which some sources are randomly removed according to 
the probability of their P(S) nd sources separated by less than 50" are merged to a single 
source. 

The simulated photons are isotropic in this area except edge part. We model the Galactic 
diffuse background using the models (gal_2yearp7v6_v0.fits) recommended by the LAT team, 
although it only contributes an uniform background to the stacked image. The EGB is not 
modeled, because the photons from point sources are actually regarded as the EGB. 

To avoid the edge effect, we only stack the sources with the coordinates between right 
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ascensions of 100° and 265°, decimations of -5° and +65°. Only about 2.2 xlO 5 of randomly 
selected simulated sources are stacked to obtain the similar flux with the stacked source of 
the FIRST sources. Other ~ 1 x 10 5 random positions are also stacked to obtain the similar 
photon number density with the stacked image of the FIRST sources. The photon number 
density profile of the stacked image is shown in Figure 9. Then we use the ML method to 
obtain the mean flux, which is 1.08 xlO -12 photon cm~ 2 s _1 . Its TS is 41.9, and \ 2 is 31.7. 
The 1 a error is 1.5xl0~ 13 photon cm~ 2 s _1 . Therefore, the derived mean flux is equal to 
the input flux within 1 a error, indicating that our flux estimation is correct and reliable. 

The stacked image contains some regions which are not independent, leading to most 
photons being counted more than once. If a photon comes from other point source rather 
than the center source, it will be regarded as a part of the backgrounds and does not affect 
our results, as proved by the simulation. 

In order to investigate the effect of the complex source merging and the sidelobe remov- 
ing on our results, we use different 6 cu t to merger the complex sources. The result, shown in 
Figure 10, presents that the fraction decreases when the 8 cut increases. This is due to two 
factors: (1) many sidelobes and complex sources are merged; (2) many independent sources 
are wrongly merged. 

For more detail, we also merger the simulated sources using various 6 cut to calculate the 
total flux of the samples, in which the Galactic diffuse background is not added. The result is 
shown in Figure 10. The slope of the first curve is obviously steeper than that of the second 
curve when 9 cut < 100", but the slopes of two curves are nearly the same when 9 cut > 100". 
This indicates that the effect of sidelobe and complex source merging is obvious only when 
@cut ;$ 100". However, when Q cut = 50", the effect of the complex source merging and the 
sidelobes removing is small, and only makes the fraction to be overestimated by about 2 %. 
The wrong merging of independent sources makes the fraction to be also underestimated by 
about 2 %. Therefore, the cutoff of 50" is appropriate and only introduces a small uncertainty. 

Consideri ng that the FIRST is a flux limited sample of radio sources and incomplete at 



the faint limit ( IWhite et al.lll997l ). we consider that when the sample is more complete or its 
radio flux limit further decreases, the contribution of the FIRST sources to the EGB should 
increases, but the exact fraction is not clear because it depends on the shape of the radio 
source count distribution below the FIRST flux limit and the correlation between the source 



radio and 7-ray lu minosity. Nevertheless , normal galaxies (Bhatt acharva fc Sreekumarll2009 



Fields et al.ll2010f ) or radio-quiet AGNs (llnoue et al.ll2008 



Inoue fc Totani 



2009), which can- 



not be well traced by the FIRST survey (only a few fraction of those sources can be included 
in FIRST), can also contribute a fraction of the EGB, we think that the point sources can 
contribute most of the EGB. 
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5. Conclusions 

In the paper, we use the method of image stacking to study the origin of the EGB, and 
find that the FIRST sources undetected by LAT can contribute about (56±6) % of the EGB. 
Considering the flux limit and incompleteness of the sample at the faint limit, we think that 
most of the EGB is distributed by point sources which are not resolved by LAT because of 
source confusion and weak flux. The main contributors of the EGB maybe blazars, non- 
blazars AGNs and starburst galaxies. But it is difficult to derive the exact fraction of each 
population contributing to the EGB using our method alone. 
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Fig. 1. — The normalized source number density profiles around the FIRST sources with 
Fpeak >100 mJy, for all sources by dot-dashed (blue) line, for the sources randomly removed 
with the probability of their P(S) by dotted (magenta) line, and for the sources with P(S)<0.1 
by dashed (green) line. The solid (red) line represents the uniform density profile. 
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Fig. 2. — The mean number distribution of FIRST sources around one FIRST source. The 
dashed (red) line represents an estimate of the chance coincidence rate. 
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Fig. 3.— The photon number density profile of the stacked FIRST sources with 
Fpeak >1 rnJy. The dashed (green) line represents the diffuse background, while the dotted 
(magenta) line represents the point source added by an arbitrary factor. The solid (red) 
line represents the best-fit model. The x 2 over degree of freedom (dof) is x 2 /dof = 30.7/31, 
indicating that the model gives a satisfactory fit to the data. The TS is 43.0, corresponding 
to a significance of ~ 6.7a. 
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Fig. 4. — o e {&) for the stacked image of the FIRST sources. 



-18- 



2.21 - 



CM 


2.209 


CD 
T3 




O 


2.208 


CM 
CD 

■o 





2.207 



2.206 - 



2.205 L 





0.2 



0.4 0.6 

6[deg] 



0.8 



Fig. 5. — The photon number density profile of the stacked imaginary sources. 
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Fig. 6. — Same as Figure 5, but angular range is extended to 2° and a small number of 
imaginary sources is stacked. In the figure, the effect of the nearby 2FGL sources is more 
obvious. 



CO 



CM 



o 



CD 

&_ 

Z3 
CO 

o 

Q. 
X 
LU 



1.15 



1.1 



1.05 



0.95 



0.9 



10' 



-20- 



Stacked sources 
Simulated sources 




10' 



10' 



Energy[GeV] 



Fig. 7. — The exposures of the stacked source and the simulated source with coordinates of 
(RA, DEC) = (190°, 30° ). 
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Fig. 8. — The photon number density profile of the simulated source with coordinates of 
(RA, DEC) = (190°, 30° ). 
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Fig. 9. — The photon number density profile of the stacked image using the simulated data. 
The dashed (green) line represents the diffuse background which is uniform. The solid (red) 
line represents the best-fit model, in which TS is 41.9, and x 2 is 31.7. 
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Fig. 10. — The percentage of the EGB resolved as a function of angular cut in the FIRST 
catalog. The dashed (red) line represent the expected dependence of the percentage on the 
angular cut for a random source distribution, which is obtained from the simulation. For 
comparison, two curves are normalized to have same value at 9 cut = 120". 



